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ABSTRACT 

The observational determination of the behaviour of the star formation rate (SFR) 
with look-back time or redshift has two main weaknesses: 1 - the large uncertainty of 
the dust/extinction corrections, and 2 - that systematic errors may be introduced by 
the fact that the SFR is estimated using different methods at different redshifts. Most 
frequently, the luminosity of the Ha emission line, that of the forbidden line [Oll]A3727 
and that of the far ultraviolet continuum (UV) are used with low, intermediate and 
high redshift galaxies respectively. 

To assess the possible systematic differences among the different SFR estima¬ 
tors and the role of dust, we have compared SFR estimates using Ha, SFR(Ha), 
[Oii]A 3727A, SFR(OH), UV, SFR(UV) and FIR, SFR(FIR) luminosities of a sample 
comprising the 31 nearby star forming galaxies having high quality photometric data 
in the UV, optical and FIR. 

We review the different “standard” methods for the estimation of the SFR and 
find that while the standard method provides good agreement between SFR(Ha) and 
SFR(FIR), both SFR(OH) and SFR(UV) are systematically higher than SFR(FIR), 
irrespective of the extinction law. 

We show that the excess in the SFR(OH) and SFR(UV) is mainly due to an 
overestimate of the extinction resulting from the effect of underlying stellar Balmer 
absorptions in the measured emission line fluxes. Taking this effect into consideration 
in the determination of the extinction brings the SFR(OH) and SFR(UV) in line with 
the SFR(FIR) and simultaneously reduces the internal scatter of the SFR estimations. 

Based on these results we have derived “unbiased” SFR expressions for the 
SFR(UV), SFR(OH) and SFR(Ha). We have used these estimators to recompute the 
SFR history of the Universe using the results of published surveys. The main results 
are that the use of the unbiased SFR estimators brings into agreement the results 
of all surveys. Particularly important is the agreement achieved for the SFR derived 
from the FIR/mm and optical/UV surveys. The “unbiased” star formation history of 
the Universe shows a steep rise in the SFR from z = 0 to z = 1 with SFR oc (I -I- 
followed by a decline for 2 > 2 where SFR oc (1 -I- z)~^'^. Galaxy formation models 
tend to have a much flatter slope from z = 0 to z = 1. 

Key words: Stars: formation. Galaxies: star forming, Galaxies: HH, Galaxies: evo¬ 
lution 


1 INTRODUCTION 

The knowledge of the history of star formation at cosmic 
scales is fundamental to the understanding of the forma¬ 
tion and evolution of galaxies. Madau and collaborators 
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combin ed th e results from the ultraviolet surveys of Lilly 
et al. (1996) with the information from the Hubble Deep 
Field to give an estimate of the star formation history from 
z = 0 to z = 4. Subsequent studies have indicated the cru¬ 
cial role played by dust in the estimates of SFR. Large cor¬ 
rection factors were suggested for z > 1-2 by several au¬ 
thors (Meurer et al. 1997, Meurer, Heckman and Calzetti 
1999, Steidel et al. 1999, Dickinson 1998). But even these 
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large dust extinction corrections do not seem to be enough 
to bring the optical/UV SFR estimates in line with the 
mm/sub-mm ones (Hughes et al. 1998; Rowan-Robinson 
et al. 1997 (RR97); Chapman et al. 2001). 

On top of the uncertainties associated with the extinc¬ 
tion correction, most of the SFR estimates have been per¬ 
formed using expressions derived from spectra constructed 
using population synthesis methods, an approach that re¬ 
quires four rather uncertain ingredients: 1) an initial mass 
function (IMF); 2) a stellar evolutionary model grid giving 
the luminosity and effective temperature as a function of 
time; 3) a stellar atmospheres grid that assigns a spectrum 
to each star for a given luminosity and effective temper¬ 
ature, and 4) a star formation history. The fact that the 
redshift evolution of the SFR is constructed using different 
estimators at different redshift ranges is a potential source 
of systematic effects with redshift that can distort the shape 
of the evolutionary curve. In practice the Ha luminosity is 
used to estimate the SFR for galaxies with redshifts up to 
0.4; the [OH] A3727A line , for those with 0.4>z>1.0 and the 
UV continuum luminosity for galaxies with z>2.0. In addi¬ 
tion, dust extinction corrections are not treated uniformly 
over the whole redshift range. 

It is therefore important to ensure that there are no 
systematic differences between the different estimators and 
corrections that can distort the results. 

This paper has two main aspects, in the hrst 5 sections 
we review the “standard” methods for the estimation of the 
SFR and test the consistency of the different SFR estimators 
by applying them to a sample of well studied nearby star 
forming galaxies and comparing the results. In the absence 
of systematic differences among them, all should give the 
same SFR for each one of the galaxies in the sample. We 
then use the results of the nearby sample to construct a set 
of “unbiased” SFR estimators and apply them to published 
surveys. 


2 THE REFERENCE GALAXY SAMPLE 


To critically test the possibility of systematic differences be¬ 
tween the SFR estimators we have compiled from the lit¬ 
erature a sample of all the nearby well studied star form¬ 
ing galaxies for which good data is available in Ha, H/3, 
[Oll]A3727, UV continuum and FIR. We call it the “refer¬ 
ence” sample of galaxies. 

The resulting sample consists of galaxies classified either 
as HII galaxy, Starburst (SB) or Blue Compact (BCG) and 
although it covers a range of galaxian properties, a large 
fraction of them is of low luminosity and low metal content. 
Work by Koo and collaborators (Koo et al. 1996, Lowenthal 
et al. 1997) has shown deep similarities between the faint 
blue galaxies found typically at z ^ 0.5 and the type of 
galaxies in our sample. This could also be the case for Lyman 
limit galaxies (Giavalisco, Steidel & Macchetto 1996). 


Spectroscopic data for 14 of the galaxies come from Mc- 
Quade, Calzetti and Kinney ( 1995| ) where they combined 
spectral data in the UV from the lUE (Kinney et al. 1993) 


with optical observations. The optical spectra were obtained 
with the KPNO 0.9m telescope covering a spectral range 
from 3500 A to 8000 A with 10 A resolution. The circular 


aperture used by McQuade et al. (1995), 13.5” in diame¬ 


ter, matches the 10” x 20” lUE aperture. Therefore no area 
renormalization was necessary for these galaxies when com¬ 
paring UV with optical data. 

Optical data for other 12 galaxie s wer e obtained by 
Storchi-B ergmann, Kinney a. nd Challis ( 1995| ) and combined 
with the (Kinney et al. 1993) lUE ultraviolet atlas. The op¬ 
tical observations were made with the Im and 1.5m CTIO 
telescopes. They used a long slit of 10 arcsecond width which 
was vignetted to equal the 10” x 20” lUE aperture. The Im 
telescope covered a range from 3200 to 6400 A with a res¬ 
olution of 5.5 A and the 1.5m telescope covered from 6400 
to 10000 A with a resolution of 8 A. In general their spec¬ 
tra flux levels agree within 20% with the spectra from the 
lUE satellite. In the cases where the differences were larger 
(observations made under non-photometric conditions) they 
assumed that the flux given by the lUE observations is the 
correct one. 

The optical spectra of CAM0840, TOL1247 and 
CAM1543 were observed by Terlevich et al. (1991). 
TOL1247 was observed under photometric conditions with 
the 3.6m ESO telescope and an 8x8 arcsec aperture. The 
spectra of CAM0840 and CAM1543 were obtained at Las 
Campanas Observatory with the 2.5m telescope. An aper¬ 
ture of 2x4 arcsec and a resolution of about 5 A were used. 
The ultraviolet spectra of these galax ies ar e from lUE satel¬ 
lite observations, Terlevich et al. (1993) (CAM0840 and 
TOL1247) and Meier & Terlevich ( |1981| ) (CAM1543). 

The optical data of MRK309 is from the UCM 
objective-prism survey (Gallego et al. 1996). The UV data 
comes from the obse rvatio ns made with lUE and reported 
by McQuade et al. (1995). 

Table ^ shows the complete sample of star forming 
galaxies used in this study. All the objects are located at 
large galactic latitudes (jbj >25). The a ppare nt blue mag¬ 
nitudes are from De Vaucouleurs et al. (1991) and are cor¬ 
rected for galactic extinction, internal extinction and aper¬ 
ture differences. Radial velocities were corrected by the Lo¬ 
cal Group motion using the NED0 velocity correction soft¬ 
ware (Table ^). 

The 60 n m IRA S data for all the galaxies come from 
Moshir et al. ( 1990 ). 

We computed the equivalent width of H/3 using the eas¬ 
ily deduced expression. 


EW{H(3) = 

Le(4861A) (1) 

= 2.5 x 10-32 io(o.4Mb) 

where L{H(5) and Lc(486lA) are the H/3 and adjacent con¬ 
tinuum luminosities respectively. The field of view used to 
estimate Mb is larger than the apertures for the spectro¬ 
scopic observations, therefore the estimated EW(H/3) rep¬ 
resents a lower limit. Where Mb was not available, as for 
CAM0840, CAM1543 and TOL1247, we used EW(H/3) from 
the spectrophotometry of Terlevich et al. (1991). 

The extinction correction to the observed fluxes both on 
the continuum and on the emission lines was estimated fol¬ 
lowing standard procedures. Two different extinction curves 
were used: the Milky Way extinction law (MW) given by 
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Name 

Type 

ms 

Velocity 
( km s“^) 

F(Ha) 

F(H/3) 

F(H7) 

F(OII) 

Fa(UV) 

(*) 

F„ 

(60 /J,m) 

Ref. 

NGC7673 

HII 

12.86 

3673.20 

608.0 

109.8 

22.1 

446.6 

14.6 

4.95 

(a) 

CAM0840 

HII 

- 

9000.00 

127.1 

37.1 

15.8 

46.1 

3.6 

0.29 

(b) 

CAM1543 

HII 

- 

11392.11 

191.0 

52.5 

22.6 

31.3 

4.1 

- 

(e) 

TOL1247 

HII 

- 

14400.00 

504.5 

134.9 

58.7 

159.2 

11.0 

0.51 

(b) 

NGC1313 

HII 

9.29 

269.86 

148.4 

17.3 

3.5 

105.4 

6.60 

14.56 

(d) 

NGC1800 

HII 

12.87 

614.07 

125.9 

26.0 

4.0 

199.3 

15.00 

0.77 

(d) 

ES0572 

HII 

14.16 

871.51 

647.9 

106.7 

52.8 

283.7 

19.95 

0.86 

(d) 

NGC7793 

HII 

9.37 

253.35 

221.9 

33.0 

- 

110.6 

9.37 

8.89 

(d) 

UGCA410 

BCDG 

15.45 

854.20 

324.4 

80.2 

21.5 

104.6 

10.4 

0.30 

(a) 

UGC9560 

BCDG 

14.81 

1305.00 

529.4 

144.6 

50.7 

334.5 

18.4 

0.71 

(a) 

NGC1510 

BCDG 

13.45 

737.95 

512.1 

116.8 

31.6 

340.9 

16.52 

0.89 

(d) 

NGC1705 

BCDG 

12.58 

401.23 

434.1 

130.0 

22.3 

265.0 

93.67 

0.87 

(d) 

NGC4194 

BCG 

12.86 

2598.46 

1946.6 

239.1 

70.1 

385.3 

14.2 

23.52 

(a) 

IG1586 

BCG 

14.74 

6045.25 

229.3 

38.8 

8.6 

138.8 

5.0 

0.96 

(a) 

MRK66 

BCG 

15.00 

6656.66 

121.3 

43.2 

11.0 

148.2 

9.9 

0.54 

(a) 

Harol5 

BCG 

- 

6498.24 

301.1 

81.0 

25.7 

264.6 

18.15 

1.35 

(d) 

NGC1140 

BCG 

12.56 

1503.13 

1400.0 

350.3 

121.3 

1010.0 

44.15 

3.36 

(d) 

NGC5253 

BCG 

10.87 

270.61 

7717.0 

2406.0 

973.2 

4370.3 

99.05 

30.51 

(d) 

MRK542 

BCG 

15.80 

7518.56 

88.9 

17.1 

- 

3.8 

5.75 

0.48 

(a) 

NGC6217 

SB 

11.66 

1599.81 

607.4 

92.9 

14.8 

108.1 

15.3 

11.05 

(a) 

NGC7714 

SB 

12.62 

2993.82 

2795.8 

539.9 

196.9 

951.1 

26.5 

10.44 

(a) 

NGC1614 

SB 

13.28 

4688.15 

1069.4 

92.0 

19.8 

95.2 

4.8 

32.71 

(a) 

NGC6052 

SB 

13.40 

4818.45 

565.0 

122.7 

39.6 

376.2 

9.9 

6.31 

(a) 

NGC5860 

SB 

14.21 

5532.07 

296.7 

26.8 

- 

15.4 

5.5 

1.64 

(a) 

NGC6090 

SB 

14.51 

8986.89 

675.3 

123.7 

42.0 

153.6 

9.6 

6.45 

(a) 

IG214 

SB 

14.16 

9161.10 

152.2 

21.4 

5.4 

35.9 

6.2 

5.22 

(a) 

MRK309 

SB 

14.61 

12918.07 

108.0 

16.2 

- 

4.4 

2.4 

3.43 

(c) 

NGC3049 

SB 

12.77 

1321.32 

513.1 

116.1 

51.9 

148.5 

10.25 

2.82 

(d) 

NGC4385 

SB 

12.90 

1981.44 

950.0 

150.7 

51.2 

261.0 

11.82 

4.73 

(d) 

NGC5236 

SB 

7.98 

304.31 

4507.0 

940.1 

154.7 

440.0 

185.47 

110.30 

(d) 

NGC7552 

SB 

11.13 

1571.17 

2064.0 

277.8 

44.5 

243.8 

19.97 

72.03 

(d) 


Table 1. Observed properties of the galaxy sample. The blue apparent ma gnitu des are from De Vauco uleurs et al. (1991). The e mission 


(b) T erlevi ch et al. (|l993|), (c) Galleg o et al. ( 1996[) 


lines and the ultravi olet c ontinuum fluxes are from (a) McQu ade et al. 
and McQuade et al. (1995), (d) Storchi-Bergmann et al. (1995), (e) Meier & Terlevich (1981) and Terlevich et al. (1991). Units for the 
intensities are 10“^® erg s“^ cm“^ except (*) where the units are 10“^® erg s“^ cm“^A“^ 

Janskys. 


The FIR fluxes are IRAS at 60 fj.m in 


Seaton ( 1971: ) and Howarth ( 1983|) and the L arge Magel¬ 
lanic Cloud one (LMC) given by Howarth ( 1983| ). A detailed 
description of the procedure can be found in Appendix 1. 


SFR(Ha){MQyr-^) = 7.9 x 10“'‘^Lrrc(erg s“^ ) (2) 

valid for a Te = 10"^K and Case B recombination, i.e. all the 
ionizing photons are processed by the gas. 


3 STANDARD SFR ESTIMATES 

In this section we describe the four estimators we have used 
to compute the SFR from the luminosities of the Hq and 
[Oll]A3727 nebular lines, the UV luminosity and the FIR 
continuum. The expressions used here are the ones we found 
most frequently referred to in the literature. All expressions 
are for a Salpeter IMF (N(m)oc with masses varying 

from 0.1 to 100 M©, solar metallicity and continuous star 
formation. To convert the fluxes into luminosities we used a 
Hubble constant of 70 km s“^ Mpc“^. 


3.1 From Ha Luminosities. 

For the estimate of the SFR from Ha luminosities we used 
t he ex pression given by Kennicutt, Tamblyn and Congdon 


3.2 From [Oll]A3727 Luminosities. 


The doublet [Oll]A3727 luminosity is used as a SFR tracer 
for objects with redshift larger than 0.4 where Ha is shifted 
outside the optical range. Unlike Ha, the [Oll]A3727 inten¬ 
sity depends not only on the electron temperature and den¬ 
sity but also on the degree of ionization and on the metallic¬ 
ity of the gas. In practice a semi-empirical approach is used 
combining the SFR(Ha) with the average Ha / [OlI]A3727 
ratio given by Gallagher et al. (198^ using a sample of 75 
blue galaxies and by Kennicutt (1992) from a sample of 90 
normal and irregular galaxies. 


SFR([Oll]){MQyr-^) = 1.4 x 10-"Lor/(erg s-^ ) (3) 
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3.3 From UV Continuum Luminosities. 

The UV continuum luminosity is used as a SFR tracer in 
objects with redshift higher than 1-2. At these redshifts all 
strong emission lines, apart from Lya, are shifted outside 
the optical range. 

In young stellar clusters, the UV spectrum is dominated 
by the continuum emission of massive stars. In evolutionary 
synthesis models of starbursts, after a short initial transient 
phase, the UV luminosity per unit frequency becomes pro¬ 
portional to the SFR, 

SFR{UV){MQyr-^) = 1.4 x 10“^®L4erg s"^ Hz~^) (4) 

This equation is valid from 1500 to 2800 A where the 
integrated spectrum is nearly flat in Ft/ for a Salpeter IMF, 
continuous mode of star formation and solar metallicity. 
This region of the UV is not affected by the Lya forest 
and the contribution from old populations is still very small 
(Kennicutt 1998, Madau et al. 1998). 

This SFR indicator though, is extremely sensitive to 
uncertainties in the reddening correction. 


4 COMPARING THE DIFFERENT SFR 
ESTIMATORS 


We have applied the commonly used SFR estimators to our 
reference sample of star forming galaxies. Given that the 
estimators are all for the same IMF and stellar models we 
do not expect these aspects to introduce any scatter. Fig¬ 
ure @ shows the SFR(Ha), SFR(OII) and SFR(UV) plotted 
against the SFR(FIR). Clearly the sample shows a correla¬ 
tion plus a large scatter. 

To simplify the analysis and simultaneously make use 
of the fact that SFR(FIR) is probably the best SFR esti¬ 
mator available, we will use in what follows FIR normalized 
SFR, i.e. the SFR relative to SFR(FIR). The FIR normal¬ 
ized SFR(Hq), SFR(OII) and SFR(UV) are: 


AHa 


= log 


SFR{Ha) 

SFR(FIR) 


AOII = log 


SFR{OII) 

SFR{FIR) 


AUV 


log 


SFR(UV) 

SFR(FIK) 


3.4 From FIR. 


To transform the observed 60/rm luminosity to SFR we as¬ 
sumed that a fraction e of the ultraviolet/optical flux emit¬ 
ted by stars is absorbed by dust and reemitted as thermal 
emission in the far infrared (10 - 300 pm). Work by, among 
others, Mas-Hesse and Kunth (1991) indicate that even for 
a small amount of reddening the fraction e is very close to 
unity. This, plus the fact that no dust extinction correction 
is necessary, justifies the assumption that the FIR luminos¬ 
ity is an excellent indicator of the total UV/optical emis- 
sivity of a galaxy. Thus, the relation between the total FIR 
luminosity and the star formation rate (Kennicutt 1998) is 

SFR{FIR){MQyr-^) = 1/e 4.5 x 10"'‘‘‘LF/H(erg s"^ ) (5) 


where the FIR lumino sity is given by Lrir ~ 1.7 x Leoum 
( Chapman et al. 200C| ). 

A difficult parameter to quantify is the fraction of ion¬ 
izing photons that escape from the nebula. Heckman et al 
(2001) show that in five of the UV brightest local starburst 
galaxies the fraction of photons escaping is less than about 
6% while Steidel, Pettini and Adelberger (2001) claim a 
higher escape fraction. Bearing in mind that the two samples 
differ, these results may not be contradictory. Tenorio-Tagle 
et al. (1999) have shown that the escape of photons from 
a starburst may be time dependent with a very large es¬ 
cape probability during the most luminous phases and little 
escape at other stages. 

It is important to point out that if the dust and ion¬ 
ized gas distributions are similar, i.e. they coexist spatially, 
the FIR and the emission line luminosities will be similarly 
affected by the escape of photons. Thus, under this condi¬ 
tion, the ratios of emission line fluxes to FIR flux are, to 
first order, independent of the fraction of escaped photons 
and therefore not very sensitive to variations in the photon 
escape from nebula to nebula. 


This is better seen in the distribution histograms of the 
normalized SFR as shown in Figure As reference we in¬ 
cluded in parts a,d,g of the figure the normalized SFR com¬ 
puted using the observed luminosities. The central and right 
columns show the dust extinction corrected ratios using the 
MW and the Calzetti extinction laws respectively. The cor¬ 
rections were applied following the common methodology 
and are described in Appendix 

Our main conclusion is that irrespective of the extinc¬ 
tion law applied, the SFR(Ha) is close to the SFR(FIR) 
while both SFR(OII) and SFR(UV) show a clear excess. 
The excess is much larger for SFR(UV) than for SFR(OII) 
suggesting a wavelength dependent effect, probably an ex¬ 
tinction over-correction. Bearing in mind that our reference 
sample has a large fraction of low metallicity and low extinc¬ 
tion galaxies this result suggests that applying these stan¬ 
dard methods to estimate SFR will systematically overesti¬ 
mate the SFR in samples at intermediate and high redshifts 
where either SFR(OII) or SFR(UV) are used. This result 
is in apparent contradiction with what has been found and 
shown in Madau-type plots in recent years, where the SFR 
obtained from UV and optical data are much lower than 
that obtained from mm and sub-mm observations at inter¬ 
mediate and high redshifts. In order to reach agreement be¬ 
tween both determinations, fixed (and somehow arbitrary) 
amounts of extinction have been applied to the UV/optical 
data, because at the moment, the intermediate and high 
redshift samples do not allow a reliable determination of the 
dust extinction. It is worth noting that Steidel et al. (1999) 
applied a fixed correction to their sample that is close to the 
average AUV. On a positive note we should indicate that 
the application of the reddening corrections reduce consid¬ 
erably the scatter in all three normalized SFR estimators as 
we will show below. 

It is important to clarify the origin of the detected ex¬ 
cess in the extinction corrected AOII and AUV. There is one 
effect that is not normally taken into account, namely that 
the presence of an underlying young stellar population with 
deep Balmer absorptions will bias the observed emission line 
ratios towards larger Balmer decrement values, mimicking 
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the dust extinction effect. In the next section we will re- 
calcnlate the SFR for the different tracers bnt including an 
estimation of the effect of an underlying population. We will 
also estimate the effect of photon escape in AUV. 


5 UNDERLYING POPULATION AND 

PHOTON ESCAPE 

5.1 Underlying absorption corrections to the 
emission line fluxes 

A clear signature of a population of young and intermedi¬ 
ate age stars is the presence of the Balmer series in ab¬ 
sorption in their optical spectrum. A complication is that 
in star forming objects, the Balmer emission lines from the 
ionized gas appear superimposed to the stellar absorption 
lines. This effect, growing in importance towards the higher 
order Balmer lines, is illustrated in Figure ^ where an ex¬ 
ample (NGC 1510) is shown. It can be seen that, while H/3 
emission is moderately affected by the absorption, all of the 
H5 emission is lost into the absorption. The equivalent width 
of the Balmer absorptions peaks at H5 - He and there is no 
detection of any absorption in Ha. This is due to the facts 
that the Ha absorption equivalent width is much smaller 
than that of H/3, and that the wings of the Ha absorption are 
difficult to detect due to the presence of forbidden [Nil] dou¬ 
blet emission at A 6548A and A 6584A, right on top of both 
wings. In spectra of poorer S/N or lower spectral resolution 
than that of Figure M, the wings of the Balmer absorptions 
are not detected and the result is an underestimate of the 
emitted fluxes and, more important for luminosity determi¬ 
nations, an overestimate of the internal extinction (Olofsson 
1995). 

The observed ratio between two emission lines (e.g. Ha 
and H/3), when the underlying absorption is included, is: 


( 6 ) 


F{Ha) F+ (Ha) - F- {Ha) 

Ffm ~ F+{Hf3)-F.{Hfi) 

where F^{Ha) and F+{HP) are the intrinsic emission line 
fluxes and F-{Ha) and F-{HP) are the intrinsic fluxes of 
the corresponding absorption lines. This expression is cor¬ 
rect in the case that the emission and the absorption lines 
have approximately equal widths. 

Including the relation between the equivalent width, the 
flux of the continuum and the intensity of the line in equa¬ 
tion ^ we obtain, 


p/rr \ 2 86 — PO^ttcJJLsl 2.86 

F{Ha) FciHp) 


_ „„ EW+{H0) 


l-Q 


(7) 


where, Fc{Ha) and Fc{Hl3) are the continuum in Ha and 
H/3 respectively, EW+ and EW_ are the equivalent widths 
in emission and in absorption respectively for the different 
lines, Q = £iv^(g[ 3 ) ratio between the equivalent 

widths of H/3 in absorption and in emission, P = [Eiv~(g/ 3 ) 
the ratio between the equivalent widths in absorption of Ha 
and H/I and P+(Ha)/ P+(H/3)=2.86 is the theoretical ratio 
between Ha and H/3 in emission for Case B recombination 
( Osterbrock 198^ ). 

The value of P can be obtained fro m spec tral evolution¬ 
ary calculations like those of Olofsson (1995). For the case 


of solar abundance and stellar masses varying between 0.1 
and 100 Mq within a Salpeter IMF, the value of P changes 
between 0.7 and 1 for ages between 1 and 15 million years 
respectively. This variation in the P parameter produces a 
change in the estimated P(Ha)/P(H/3) ratio of less than 2%, 
so in what follows we asssume P=l. 

For an instantaneous burst, the ratio varies 


betwe en 0.14 and 0.26 (Mayya (1995), Leitherer & Heckman 

(H))- 


The corresponding equation for Hy and H/3 is: 
P(P7) _ 0.47 - GQ 

F\m ~ l-Q 


( 8 ) 


where G = 1^® ratio between the equivalent 

widths in absorption of Hy and H/3 and we assume for the 
respective emis sions an intrinsic ratio of 0.47 for Case B 
recombination ( [Osterbrock 1989 ). 

The evolution of the equivalent width of the Balmer 
absorption lines has bee n ana lyzed by Conzalez Delgado, 
Leitherer and Heckman (1999). In their models the param¬ 


eter G is almost constant in time and independent of the 
adopted star formation history. We fixed the value of G to 
1 as suggested by their results. 

The effect of the underlying stellar absorptions is shown 
as a vector Q in Figure ^ (from equations ^ and |^. The 
whole time dependence is shown by the three closely grouped 
vectors. Its range is much smaller than typical observational 
errors. Dust extinction is also represented by a vector in the 
same plane (equation AS). It is possible, as these two vectors 
are not parallel, to solve simultaneously for underlying ab¬ 
sorption {Q) and extinction (Av) for every object for which 
F{H'y), F{HP) and F{Ha) are measured. 

We further illustrate the presence of underlying Balmer 
absorption in star forming galaxies in Figure H where 
we have plotted the galaxies from our sample in the log 
(F(H7)/F(H/3)) vs. log (F(Ha)/F(H/3)) plane. Also shown 
are the vectors depicting dust extinction and the underly¬ 
ing absorption. Clearly, most observational points occupy 
the region below the reddening vector and to the right of 
the Balmer absorption vector. In the absence of underlying 
absorption all points should be distributed along the extinc¬ 
tion vector. The fact that there is a clear spread below the 
extinction vector gives support to the underlying absorption 
scenario. We also find that the objects with smaller equiv¬ 
alent width of H/3 are systematically further away from the 
pure extinction vector. 

Four galaxies (NCC 3049, ESO 572, MRK 66 and 
NCC 1705) fall outside the space defined by the extinction 
and underlying absorption vectors, although two of them are 
within the errors. The other two (ESO 572 and MRK 66) 
are faint and reported as having been observed in less than 
optimum conditions in the original observations paper. 

We have used this method to estimate simultaneously 
the “real” visual extinction Av* and the underlying Balmer 
absorption Q The values of Av* were then applied to the 
UV continuum and the emission line fluxes; the corrected 
values are listed in Table 0. 


§ See Appendix ^ for a detailed discussion on the dust extinction 
corrections to the observed fluxes. 
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Figure 1. Standard SFR estimators vs. SFR(FIR). No extinction corrections were applied to the data. The solid line represents equal 
values. 


No extinction MW extinction Calzetti's extinction 



Figure 2. Histograms of the SFR rates given by the different tracers normalized to the SFR(FIR). In the left panels (a,d,g) no corrections 
were applied to the data. In the central panels (b,e,h) we corrected the Ha, [OHJA3727 and UV continuum by using the MW extinction 
curve. In the right panels (c,f,i) the Ha, [OH]A3727 and UV continuum were corrected with Calzetti’s extinction law. The median and 
standard deviation are given for each case. The number of objects is 29. 
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Figure 3. A blue spectrum of the star forming galaxy NGC 1510 is shown to illustrate the effect of stellar Balmer absorptions in the 
measurement of the emission line strengths. 


NGC1510 


NGC1614 




log[F(H«)/F(Hf?)] log[F(Ha)/F(H/?)] 


Figure 4. Logarithmic ratio of F(Hq:)/F(H/3) vs. F(H7)/F(H/3). The “observed vector” (from the intrinsic values given bv recombination 
theory to the observed ratio^can be decomposed in two vectors, one is due to pure extinction and is given by Equation |A8[ the other one 
is given by Equations^ andg and shows the effect of an underlying stellar population which is characterized by Q (see text). Plotted are 
the cases for NGC 1614 and NGC 1510. NGC 1614 is an example where the underlying absorption correction to the extinction estimate 
is not very large. NGC 1510 on the other hand, has a large correction. 
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log[F(Ha)/F(H/?)] 

Figure 5. The galaxies from our sample are plotted in the Balmer decrement plane log (F(Ho)/F(H/3)) vs. log (F(H7)/F(H/3)). The 
vectors indicate the direction of shifts due to extinction or underlying absorption from the intrinsic values given by recombination. The 
cross indicates typical Icr errors in the line ratios. Most points are, within the errors, below the reddening line and to the right of the 
underlying Balmer absorption line. 


Figure shows the result of taking into account the cor¬ 
rections for Balmer absorptions due to an underlying stellar 
population. The medians of both AOII and AHa are close 
to zero indicating that including the underlying absorption 
correction brings into agreement the SFR in the optical with 
those in the FIR. At the same time AUV shows still a pos¬ 
itive value indicating an excess with respect to the FIR es¬ 
timate. We must remember that while the ratio of emission 
line fluxes to FIR flux is not very sensitive to changes in the 
photon escape from object to object, this is not the case for 
the ratio of UV continuum to FIR fluxes. The reason being 
that while in the UV continuum we are detecting directly 
the escaped photons, i.e. those that do not heat the dust or 
ionize the gas, the emission lines and FIR fluxes are repro¬ 
cessed radiation, i.e. the product of the radiation that does 
not escape the region. 

A striking aspect is the large reduction in the r.m.s. 
scatter in the AUV from 0.70 before corrections to 0.39 af¬ 
ter corrections, i.e. about half the original value. This simple 
fact suggests the goodness of the corrections applied to the 
data. This aspect is also illustrated in Figure M when com¬ 
pared to Figure 


5.2 Photon Escape 


As discussed above, while AOII and AHa are not very sen¬ 
sitive to the fraction of escaped photons, AUV is. The es- 
cape of UV photons can be quantifie d by the parameter e 
(Rowan-Robinson & Efstathiou 1993) by assuming that 


Lfir = (FboI 
Luv = (1 — FjFsoi 


(9) 


where (1 — e) represents the escaped fraction. 

Thus, the reprocessed fraction e can be estimated 
through the observed GO/rm and ultraviolet luminosities as 


l + _kiv_ 


( 10 ) 


where the luminosities are given by 


Luvierg s"!) = Fx{UV) x 2000A 

Leo^mierg s“^) = F^leOgm) 5 10^‘^Flz 10“^^ 


( 11 ) 


The values of the UV flux in units of erg s“^cm“^A“^ 
and the fluxes at 60/rm in Jy are given in Table and the 
calculated values of the escape fraction 1 — e, in Table 

The results of including the photon escape are also 
shown in the AUV histograms of Figure |^. As discussed 
above, under the simple assumption that dust and ionized 
gas have the same spatial distribution, the values of AHa 
and AOII do not change with respect to Figure |^. 

In the left pa nel, t he UV continuum was corrected us¬ 
ing Steidel et al. ( 199£ ) simple approach of multiplying the 
observed UV flux by a fixed amount (x 5) which corresponds 
to the av erage correction fou nd in a sample of local starburst 
galaxies (Calzetti et al. 1991) . The fact that the average cor¬ 
rection for Calzetti’s sample and the average correction for 
our sample are almost identical suggests that both samples 
have been drawn from the same family of objects, i.e. they 
are similar samples. 

The central panel shows the result of using the MW 
extinction law and applying to the observed values the ex¬ 
pression A.5 where the visual extinction is Av* . The panel 
at the right shows the r esult of using the extinction law 
given by Calzetti et al. (1991). Also in this case the visual 
extinction is Av*. 

Inspection of the AUV distribution in figures |^, ^ and ^ 
shows that at least for the galaxies in the “reference” sample, 
the escape of photons is a minor effect compared to the 
underlying Balmer absorption. 
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MW extinction Calzetti's extinction 



Figure 6. The figure shows the histograms of the normal ized SFR after the underlying stellar absorption effect is deducted from the 
extinction estimates, i..e. Ay is used instead of Ay (Section 5.1). Panels a, c, e show the distributions after correcting the Ha, [OIIJA3727 
and UV continuum using the MW extinction curve, b, d, and f show the result of using Calzetti’s extinction curve. The median and 
standard deviation are given for each case. The total number of objects is 25. 



0.01 0.10 1.00 10.00 100.00 

SFR(FIR) 


Figure 7. Corrected SFR es 
photon escape (see Section 5.2 


timators vs. SFR(FIR). The corrections include the underlying Balmer absorption (see Section |5.l| ) 
The solid line represents equal values. 


and 


Comparing the MW and Calzetti’s extinction cor¬ 
rections we can conclude that the corrections using the 
Calzetti’s extinction curve give a SFR(UV) identical inside 
the errors to the SFR(FIR) plus a substantial reduction in 
the dispersion of the SFR estimates. These two points jus¬ 


tify the use of Calzetti’s extinction corrections in samples 
of starburst galaxies similar to the ones used in this work. 
It is a remarkable result that the application of our method 
succeeds in cutting down the scatter present in the origi¬ 
nal data to less than a half. These two results, i.e. the very 
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Name 

Av 

Av* 

Qx 100 


(1-e) 

EW(H/3) 

NGC7673 

1.84 

1.11 

27 

-1.50 

0.11 

4.69 

CAM0840 

0.50 

0.43 

3 

-1.26t 

0.33 

121* 

CAM1543 

0.67 

0.67 

0 

-0.70t 

- 

224* 

TOL1247 

0.75 

0.75 

0 

-0.47t 

0.46 

97* 

NGC1313 

3.06 

2.45 

23 

-0.601 

0.02 

0.03 

NGC1800 

1.47 

0.56 

33 

-1.65 

0.56 

1.10 

ES0572 

2.10 

1.87 

0 

-1.96t 

0.48 

14.81 

NGC7793 

2.38 

- 

- 

-1.34 

0.04 

0.06 

UGCA410 

0.97 

0.37 

23 

-1.84 

0.58 

36.52 

UGC9560 

0.69 

0.34 

14 

-2.02 

0.51 

36.63 

NGC1510 

1.19 

0.62 

22 

-1.71 

0.43 

8.42 

NGC1705 

0.43 

0.0 

31 

-2.42 

0.81 

4.20 

NGC4194 

2.91 

2.73 

8 

-0.26 

0.02 

10.14 

IG1586 

2.02 

1.37 

25 

-0.91 

0.17 

9.51 

MRK66 

0.00 

0.00 

23 

-1.94 

0.42 

13.50 

Haro 15 

0.73 

0.27 

18 

-1.48 

0.35 

- 

NGC1140 

0.93 

0.62 

13 

-1.78 

0.34 

11.19 

NGC5253 

0.32 

0.12 

9 

-1.33 

0.11 

16.07 

MRK542 

1.66 

- 

- 

-1.32 

0.32 

11.23 

NGC6217 

2.30 

1.46 

31 

-0.74 

0.05 

1.30 

NGC7714 

1.65 

1.57 

4 

-1.23 

0.09 

18.40 

NGC1614 

3.90 

3.47 

17 

-0.76 

0.01 

5.82 

NGC6052 

1.33 

0.98 

14 

-0.72 

0.06 

8.68 

NGC5860 

3.77 

- 

- 

-0.91 

0.12 

4.02 

NGC6090 

1.80 

1.62 

8 

-0.781 

0.06 

25.00 

IG214 

2.54 

2.08 

18 

-0.61 

0.05 

3.14 

MRK309 

2.36 

- 

- 

2.081 

0.03 

3.68 

NGC3049 

1.21 

1.01 

0 

-1.14 

0.13 

4.49 

NGC4385 

2.20 

2.10 

4 

-1.02 

0.09 

6.60 

NGC5236 

1.44 

0.56 

32 

-0.83 

0.06 

0.44 

NGC7552 

2.66 

1.85 

30 

0.48 

0.01 

2.38 


Table 2. Extinction properties, ultraviolet slope and photon 
escape probability of the galaxy sample. The second column is 
the extinction in magnitudes obtained with the MW extinction 
curve (see Appendix) without taking into account the effect of 
the underlying population. The third column is the extinction 
in magnitudes obtained by combining the MW extinction curve 
with the effect of an underlying population. The ratio between 
the equivalent widths of H/3 in absorption and in emission as 
is given in the fourth column. Values of the 


defined in Section 5.1 


ultraviolet slope /3 (see appendix) are from Meurer et al. (1999) 
except for those ma.r ked w ith a | where the values were extracted 
from Calzetti et al. (1994) and those marked with a f where we 
calculated ourselves the paramet er 3 from the lUE data base 
spectra, following Calzetti et al. (1994) prescription. The escape 
probability e is given in the sixth column (see text). Equivalent 
widths were calculated using Equation except when marked 
with a * where they were taken from published spectroscopy. 


close agreement between the SFR(UV) and SFR(FIR) plus 
the rather small scatter in the ratio of these two estimators 
gives confidence in the use of our approach to estimate dust 
extinction and star formation rates in starburst galaxies. 


6 UNBIASED SFR ESTIMATORS 


The central result we have found is that the extinction cor¬ 
rection including the effects of an underlying stellar Balmer 
absorption brings into agreement all four SFR estimators, 
and that the photon escape correction seems to play a mi¬ 
nor role. 

We have seen that the SFR given by the UV continuum 
is equal to the SFR given by the FIR if the extinction cor¬ 
rection is estimated using the Calzetti extinction curve and 
including the underlying Balmer absorption corrections. It 
is reassuring that the simple inclusion of the underlying ab¬ 
sorption correction brings into agreement theory and obser¬ 
vations. We have used the results from the previous sections 
to obtain SFR estimators that are free from these systematic 
effects. 


SFR{Ha){MQyr-^) 

= 1.1 X 10"'‘^ 

Life, (erg s“^ ) 

SFR{[Oll]){Meyr-^) 

= 8.4 X 10"'‘^ 

Lon (erg s“^ ) 

SFR{UV){MQyr-^) 

= 6.4 X 10"^® 

L^{eTg s“^ Hz 


where the SFR(Hq:) is obtained from the reddening cor¬ 
rected Ha luminosity while the expressions for the SFR(OII) 
and the SFR(UV) are for the observed fluxes. 

These estimators, (the SFR(OII) and the SFR(UV) 
from equation should be applied to samples where it is 
not possible to determine the extinction. If the objects have 
similar properties to our selected sample, then any system¬ 
atic difference between the different estimators should be 
small. 

We have applied this new set of calibrations to the SFR 
values given by different authors at different redshifts: Gal- 
lego et al. (1995) Ha survey for the local Universe, Cowie 
et al. (199E) [0II]A3727 sample and Lilly et al. (1995) UV 
continuum on e bet ween redshifts of 0.2 and 1.5 and Con¬ 
nolly et al. ( 19971) be tween 0.5 and 2; and the UV points 
by Madau et al. (1996) and Steidel et al. ( 199£| ) at redshift 
higher than 2.5. We give in Table ^ the complete list of sur¬ 
veys of star formation at different redshifts that we have 
used plus their different tracers and computed SFR. 

The results are plotted in Figure ^ We have also plot¬ 
ted the results obtained by Hughes et al. (1998) based on 
SCUBA observations of the Hubble Deep Field, those by 
Chapman et al. (2001) based on sub-millimeter observations 
of bright radio sources and those by Rowan-Robinson et 
al. (1997) based on observations at 60 fim of the Hubble 
Deep Field. Clearly the mm/sub-mm results and our “unbi¬ 
ased” results agree within the errors. The “unbiased” history 
of star formation is characterized by a large increase of the 
SFR density from z~0 to z~l (a factor of about 20) and a 
slow decay from z~2 to z~5. 

The final corrected values are similar to other published 
results (e.g. Somerville, Primack and Faber 2001). But the 
fact that our procedure removes systematic (or zero point) 
differences between the different estimators, implies that the 
shape of the curve, and therefore the slopes between z = 
0 and 1 and « = 1 and 4 are now better determined. 
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X 5 MW extinction Calzetti's extinction 



Figure 8. Histograms of the SFR given by the UV continuum normalized to the SFR(FIR). In the left panel we show the original 
data with the fixed correction given by Steidel et al. (1999). The cent ral and right panels show the AUV corrected for dust extinction 
computed including the underlying Balmer absorption (see Section 5.1) and photon escape. The median and standard deviation are given 
for each case. The total number of objects is 25. 
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Figure 9. The SFR density as a function of redshift. The solid 
and open symbols represent values corrected and uncorrected for 
reddening respectively, except for the Hq values (the two lowest 
z points) which are all reddening corrected. The SFR in all cases 
are estimated using equations g, ^ Q and for open symbols 
and equations for the filled symbols. The open circle is the 
SFR based on Ho; emission line corrected by extinction from Gal¬ 
lego et al. (1995) sample. The open stars are the S FR de termined 
by [OIIJA3727 from the sample of Cowie et al. (1995) without 
any co rrectio n. The UV continuum sa mples are from Connolly 
et al. ( 1997) ftr iangles) Madau et al. ( 1996[ ) (squares) and Stei¬ 
del et al. (1999) (upside down triangles ). Th e upward arrow is 
the lower limit given by Hughes et al. (1998) based on sub-mm 
observations of the Hubble Deep Field. The crosses at z = 0.55 
and 0.85 correspond to the SFR estimates by Rowan-Robinson 
et al. (1997). The dashed line is the lower limit given by Chap¬ 
man et al. (2001) 


7 THE HISTORY OF STAR FORMATION OF 
THE UNIVERSE 

Let’s now analyze with some detail the behaviour of the dust 
corrected curve of Figure]^. 

Figure ^ shows the corrected curve together with some 
representative theoretical models; (the details of the differ- 

© 2001 RAS, MNRAS 000, EHi 


ent sets of data are listed in table Our basic assumption 
is that the main results from our sample of galaxies in the 
nearby universe are also applicable to samples of starforming 
galaxies at intermediate and high redshift. 

An important result is that within observational errors, 
our calibrations succeed in bringing into reasonable agree¬ 
ment the FIR determinations with the optical/UV ones. A 
first implication of this agreement is that large dust extinc¬ 
tion corrections are not favoured. It can be seen that while 
the general trend present in previous determinations of the 
star formation history is preserved, the rise from z = 0 to 
2 ; = 1 is steeper than in most previous work. The procedure 
we have used minimizes any systematic differences between 
the different redshift ranges and implies that, provided the 
composition of the sample of objects at intermediate and 
high redshift is similar to that of our nearby sample, the 
shape of the curve would be preserved. 

The general behaviour of the star formation history of 
the universe can be represented with a steep rise in the SFR 
rate from 2 = 0 to 2 = 1 with SFR oc (1 + zA'^ and a decline 
from 2 = 1 to 2 = 4 where SFR oc (l-|- 2 )“^ ® (see solid lines 
in Figure |^. The steep increase to 2 = 1 is similar to 
the slope found by Glazebrook et al. (1999) using the Ho? 
luminosities only. 

We are conhdent that the procedure presented in this 
paper guarantees the removal of systematic differences be¬ 
tween the star formation estimators when they are applied 
to samples qualitatively similar to ours. This in turn means 
that the general shape and in particular the slope from 2 = 0 
to 2 = 1 and from 2 = 1 to 2 = 4 are well determined, al¬ 
though the exact value of the 2 = 0 to 2 = 1 slope is still 
very dependent on the exact value of the SFR in the nearby 
universe. 

We also included in Figure the results from some 
representative semi-analytical models of galaxy formation 
corresponding to the hierarchical cold dark matter scenario. 
The dotted lines correspond to Somerville et al. (2001) 
rCDM models with two different normalizations while the 
dashed line corresponds to model C from Baugh et al. 
(1998). 

Whatever the local value of the SFR is, the compari¬ 
son with these models of galaxy formation shows a general 
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z 


Figure 10. The unbiased SFR density as a function of redshift. 
The dotted lines represent the two normalizations in Somerville 
et al. (2001) rCDM models, and the dashed line is model C from 
Baugh et al. (1998). The solid lines represent SFR tx (1 + z)^ ® 
from z = 0 to 2 = 1 and SFR oc (1 + z)“l ® for z > 1 . 

agreement in the shape plus what we can call a local excess 
in the models, i.e. the models have a shallower slope between 
z = 0 to 2 = 1 than the data snggest or eqnivalently, that 
present semi-analytical models of galaxy formation seem un¬ 
able to reproduce the sharp increase in the SFR from z = 0 
to z = 1. 


8 CONCLUSIONS 

One aspect that has generated many discussions regarding 
the SF history diagram in Cosmology, is the lack of con¬ 
fidence on the reddening corrections. This is compounded 
with the low level of agreement found until now between 
the optical/UV and FIR determinations of the SFR. 

In the first part of this paper, we have investigated the 
possible systematic differences between SFR estimators by 
applying them to a sample of nearby star forming galaxies 
with good photometric data from the UV to the FIR. 

We found that the main source of systematic differences 
among the SFR rate estimators is related to the presence of 
stellar Balmer absorption in the spectrum of emission line 
galaxies. The main effect of the Balmer absorptions is to 
produce an overestimate of the reddennnig when their ef¬ 
fect is not included. We showed that taking into account 
the underlying Balmer absorptions effect in the estimates of 
reddening, removes most of the systematic differences be¬ 
tween the SFR estimators in the optical/UV and FIR. Fur¬ 
thermore the scatter of the SFR estimations is considerably 
reduced by the application of the corrections. We also found 
that the escape of photons plays a minor role compared to 
that of the Balmer absorptions. These results give renewed 
confidence to the estimates of SFR for star forming galaxies 
in general and for samples similar to the one presented here 
in particular. 

Thus, our central result is that the extinction 
correction including the effects of an nnderlying stel¬ 


lar Balmer absorption brings into agreement all four 
SFR estimators, and that the photon escape correc¬ 
tion seems to play a minor role. 

In the second part of the paper we used the average 
results for our sample to construct a set of “unbiased” SFR 
estimators. These “unbiased” SFR estimators expressions 
include statistically the underlying Balmer absorption and 
photon escape corrections to the extinction estimates and 
bring the four SFR estimators studied here into the same 
system. We thus obtained consistent results between the 
SFR estimators in the optical/UV and FIR. The applica¬ 
tion of these “unbiased” SFR estimators to a compilation of 
surveys has produced a SFR history of the universe where all 
surveys results agree whitin the errors. Particularly impor¬ 
tant is the level of agreement achieved between the FIR/mm 
and optical/UV SFR results. 

Our “new” and unbiased SFR history of the universe 
shows a steep rise in the SFR rate from z = 0 to z = 1 with 
SFR oc (1 -|- followed by a mild decline for z > 2 where 
SFR oc (1 -I- z)~^'^. The steep increase to z = 1 seems in 
line with recent determinations of the SFR using only the 
Ha estimator. Most galaxy formation models tend to have 
a much flatter slope from z = 0 to z = 1. 
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Redshift 

Range 

SFR 

Tracer 

log (C) 

log(SFR) 

(Standard) 

log(SFR) 

(Unbiased) 

Log Error 


References 


0.0-0.045 

Ha 

39.09 

-2.01 

-1.87 

0.2 

( 

Glalleeo et al. 1995|l 

0.1-0.3 

Ha 

39.47 

-1.63 

-1.50 

0.04 

( 

Fresse & Maddox 199^1 

0.75-1.0 

Ha 

40.22 

-0.75 

-0.88 

0.17 

( 

Glazebrook et al. 

1999^ 

0.25-0.50 

FIR 

41.85 

-1.50 

-1.50 

0.26 

( 

Flores et al. 1991 


0.50-0.75 

FIR 

42.17 

-1.18 

-1.18 

0.26 

( 

Flores et al. 1991 


0.75-1.00 

FIR 

42.49 

-0.86 

-0.86 

0.26 

( 

Flores et al. 1991 


0.25-0.50 

OH 

38.83 

-2.02 

-1.24 

0.1 

( 

Fowie et al. 199.'; 


0.50-0.75 

OH 

39.12 

-1.73 

-0.95 

0.1 

{ 

Fowie et al. 199.'; 


0.75-1.00 

OH 

39.39 

-1.46 

-0.68 

0.1 

( 

Gowie et al. 199,'^ 


1-1.25 

OH 

39.25 

-1.60 

-0.82 

t 

( 

Fowie et al. 1995 


1.25-1.50 

OH 

39.21 

-1.64 

-0.86 

t 

( 

Gowie et al. 


0.25-0.50 

UV 

25.89 

-1.96 

-1.30 

0.07 

{ 

rjllv et al. 1996 

) 


0.50-0.75 

UV 

26.21 

-1.64 

-0.98 

0.08 

{ 

rjllv et al. 1996 

) 


0.75-1.00 

UV 

26.53 

-1.32 

-0.66 

0.15 

( 

Lilly et al. 1996 

) 


0.4-0.7 

0.7-1.0 

FIR 

FIR 

42.64 

43.23 

-0.76 

-0.13 

-0.76 

-0.13 

{-H0.1)(-0.2) 

{-H0.4)(-0.2) 

I 

F 

1R97 

1R97 



0.5-1.0 

UV 

26.52 

-1.33 

-0.63 

0.15 

( 

Gonnollv et al. 

1997 

) 

1.0-1.5 

UV 

26.69 

-1.16 

-0.50 

0.15 

( 

Gonnollv et al. 

1997 

) 

1.5-2.0 

UV 

26.59 

-1.26 

-0.60 

0.15 

( 

Gonnollv et al. 

1997 

) 

2.0-4.0 

FIR 

42.80 

-0.85 

-0.85 

t 

( 

FTiip-hes et al. 199 


2.0-3.5 

UV 

26.42 

-1.43 

-0.77 

0.15 

( 

Madau et al. 199f 


3.5-4.5 

UV 

26.02 

-1.83 

-1.17 

0.2 

( 

Madau et al. 199f 


2.8-3.3 

UV 

26.28 

-1.57 

-0.91 

0.07 

( 

Dteidel et al. 1999 


3.9-4.5 

UV 

26.19 

-1.66 

-1.00 

0.1 

( 

Dteidel et al. 1999 



Table 3. SFR for surveys at different redshifts. The mean luminosity density (£) is given in units of erg s“^Mpc“^ for the case of 
the Ha and Oil emission lines and the FIR. The UV is gRen in erg s“^Hz~^Mpc“^. The SFR referred to as standard are given by the 
expressions by Kennicutt (1998) and described in Sectiori^. The SFR referred to as unbiased are the ones obtained by using Equation [l^ 
In both cases, the SFR are given in units of MQyr~^M.pc~^. 
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APPENDIX A: DUST EXTINCTION 
CORRECTIONS TO THE OBSERVED FLUXES 


Two different extinction curves are us ed: t he Milky Way 
extinc tion law (MW) given by Seaton (1979) and Howarth 
( 1983| ) and the Large Magellanic Cloud one (LMC) given 
by Howarth ( |l983|) . The main parameters of both laws are 
given in Table |Al[ 


A1 Dust Extinction Corrections to the 
Continuum Fluxes 

Calzetti and collaborators developed an empirical method to 
estimate the UV extinction (Calzetti, Kinney and Storchi- 
Bergmann, 1994). They found that the power-law index (3 
in the ultraviolet defined as Fx oc is well correlated with 
the difference in optical depth between Ha and H/3 defined 
as Tg = In ( ) where F{Ha) and F{F[I3) are the 

intensities of the Ha and H/3 emission lines respectively. This 
correlation, which is linear and independent of the adopted 
extinction law is given by, 


/3 = (1.76 ± 0.25)ts - (1.71 ± 0.12) 


(Al) 


The parameter /3 is obtained by fitting the power law to the 
lUE ultraviolet spectra. Calzetti et al. (1994) and Meurer, 
Heckman and Calzetti (1999) values of /3 are presented in 
Table as well as our estimate for CAM0840, CAM 1543, 
TOL1247, ES0572 and MRK309. 

The effect of reddening using differen t du st spatial dis¬ 
tributions can be estimated from Equation Al by comp aring 


ultraviolet with optical spectra. Calzetti et al. (1994) esti¬ 


mate the optical depth rx by solving the transfer equation 
for five different geometries, uniform or clumpy dust screen, 
uniform or clumpy scattering slab and internal dust. The 
uniform dust screen constitutes the easiest case where the 
optical depth is related to the visual extinction by 

ta = 0.921fc(A)E(R - V) 

where E(B-V) is the colour excess and fc(A) is the extinction 
law. For the other geometries, apart from the assumed ex¬ 
tinction law, the optical depth is a function of dust parame¬ 
ters such as the albedo, the phase parameter or the number 
of clumps. After comparing synthetic extinction corrected 
spectra with observations of emission line galaxies Calzetti 
et al. (1994) conclude that none of the adopted geometries 


combined with the standard MW and LMC extinction laws 
could explain the observed tight relation between Tg and 
/3 and proposed an empirical extinction law obtained from 
lUE spectra of a sa mple o f nearby starburst galaxies. 

Calzetti et al. (1994) created 6 different templates av¬ 
eraging galaxies with the same amount of dust (judging by 
their Balmer decrements). The template with Tg = 0.05 is 
taken as the reference one (free of dust). An optical depth, 
r„(A), is calculated for each template by comparing the ob¬ 
served fluxes, Fn{X} and Ei(A), 


'r„(A) = - In 


E„(A) 
El (A) 


(A2) 


where the subindex 1 corresponds to the dusty free template 
and the subindex n corresponds to the n template. For each 
template a rescaled optical depth can be defined 


© 2001 RAS, MNRAS 000, in 














An Empirical Calibration of Star Formation Rate Estimators 15 



Figure Al. Adopted extinction cu rves. The solid line is the 
empirical relation given by Cal zetti ( 1999| ). The da shed line is 
the curve for the MW (Seaton (1979) and Ho warth (1983)) and 


the dot-dashed line is the curve for the LMC (Howarth 1983). 


The ratio between the intensity of a given line F{X) and 
the intensity of H/3, F(H/3) can be expressed by: 

F{\) FRX) 


F{HP) FRHP) 


10 


-0AAv[k(X)-k(HI3)]/R 


(A7) 


where the difference fc(A) — k{Hj3) is tabulated for different 
extinction laws (Table hjl)- The total visual extinction Av, 
depends on the observed object (see Table The subindex 
o indicates the unreddened values. We use as reference the 
theoretical ratio for Case B recomb ination K)(Hq:)/J 1i(H/3~) 
=2.86 and T;(H7)/f;(H/l)=0.47 (|Osterbrock 1989|) . The 
observed flux ratios can be expressed as a function of the 
theoretical ratios and the visual extinction, 


F{HI3) 


= log 2.86 - QA[k{Ha) - k{H(5)]KY/R 
log 0.47 - QA[k{H-y) - fc(J7/?)]Av/i? 


(A8) 


This Equation was used to analyze t he p resence of an un¬ 
derlying stellar population in Section p.l| . 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science fAtjjfX style file. 


Qn(A) = 


r„(A) 


(A3) 


Averaging this quantity, Calzetti et al. (1994) found an ex¬ 
tinction curve, Q{\) which can be transformed to fc(A) (e.g. 
Seaton, 1979) by, 

fc(A) 


Q(A) = 


(A4) 


k{HI3) - k{Ha) 

where the difference k(HP) — k{Ha) is given by the Seaton 
(1979) extinction curve. The observed ultraviolet flux is re¬ 
lated to the emitted one by, 

Fobs{\) = To (A) 10-“-^"^'' (A5) 


where fc(A) is given by (Calzetti 1999) 

k{X) = -2.156 -b 1.509/A - 0.198/A^ -b 0.011/A® (A6) 

valid for the range 0.12 fim < A < 0.63 fim. 

The obtained extinction curve can be considered as 
an average of the d ifferent dust distributions described by 
Calzetti et al. (1994). 


In order to correct the ultraviolet flux using this proce¬ 
dure it is necessary to estimate Av f rom the observed Hq/H/ 3 
ratios and then apply Equation A6 to the observed ultravio¬ 


let fluxes. Physically this correction is understood assuming 
that the ionized gas is more affected by extinction than the 
stars w hich are producing the observed UV flux (Calzetti 
et al. 1994). 


No corrections were applied to the IR data. 

The different extinction curves are plotted in Figure Al 


A2 Dust extinction corrections to the Emission 
Line Fluxes 

Extinction affects the emission lines in different degrees de¬ 
pending on wavelength. Corrections are usually obtained 
from the observed ratio of Balmer lines, the intrinsic ratio, 
and an adopted interstellar extinction curve. 
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k{Ha)-k(Kl3) 

fc(H7) - fc(H/3) 

_R=fc(5464 A) 

fc([OII]A3727) 

A:(1700 A) 

MW 

- 1.25 

0.45 

3.2 

4.67 

7.80 

LMC 

- 1.18 

0.48 

3.2 

4.86 

9.54 

Calzetti 

- 0.58 

0.23 

2.7 

3.46 

5.10 


Table Al. Adopted values for the extinction curves. 
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